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The influence of tlie spin-orbit interactions on the energy spectrum of two-electron laterally cou- 
pled quantum dots is investigated. The effective Hamiltonian for a spin qubit pair proposed in 
F.BarufTa et al., Phys. Rev. Lett. 104, 126401 (2010) is confronted with exact numerical results in 
single and double quantum dots in zero and finite magnetic field. The anisotropic exchange Hamil- 
tonian is found quantitatively reliable in double dots in general. There are two findings of particular 
practical importance: i) The model stays valid even for maximal possible interdot coupling (a single 
dot), due to the absence of a coupling to the nearest excited level, a fact following from the dot sym- 
metry, ii) In a weak coupling regime, the Heitler-London approximation gives quantitatively correct 
anisotropic exchange parameters even in a finite magnetic field, although this method is known to 
fail for the isotropic exchange. The small discrepancy between the analytical model (which employs 
the linear Dresselhaus and Bychkov-Rashba spin-orbit terms) and the numerical data for GaAs 
quantum dots is found to be mostly due to the cubic Dresselhaus term. 
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I. INTRODUCTION 

The lowest singlet and triplet states of a two electron 
system are split by the exchange energy. This is a di- 
rect consequence of the Pauli exclusion principle and 
the Coulomb interaction. As a result, a spin struc- 
ture may appear even without explicit spin dependent 
interactions fi 

In quantum dot spin qubits^ the exchange interaction 
implements a fundamental two qubit gate*^"^ Compared 
to single qubit gates^i^ the exchange-based gates are 
much faster— and easier to control locally, motivating 
the solely exchange-based quantum computation.^ The 
control is based on the exponential sensitivity of the ex- 
change energy on the inter-particle distance. Manipula- 
tion then can proceed, for example, by shifting the sin- 
gle particle states electricall y^'^' or compressing them 
magnetically. ^ ^ 

The practical manipulation schemes require quantita- 
tive knowledge of the exchange energy. The configuration 
interactioufi^r— a numerically exact treatment, serves 
as the benchmark for usually adopted approximations. 
The simplest one is the Heitler-London ansatz, in which 
one particle in the orbital ground state per dot is con- 
sidered. The exchange asymptotic in this model differs 
from the exac t^'''^^ and the method fails completely in fi- 
nite magnetic fields. Extensions of the single particle ba- 
sis include the Hund-Mullikan^ Molecular Orbitalfi^iiS 
or Variational method i^^'^" Other approaches, such as 
the Hartree-FockfSi"— random phase approximation^"* 
and (spin-)density functional theory^ were also exam- 
ined. None of them, however, is reliable in all important 
regimes^ ' "'^^'^^ which include weak/strong interdot cou- 
plings, zero/finite magnetic field and symmetric/biased 
dot. 

The spin-orbit interaction, a non-magnetic spintron- 
ics workhorse?^ is a generic feature in semiconductor 



quantum dots.^^ Although it is usually weak, it may 
turn out of major importance as, for example, for the 
spin relaxatioU )^^'^^" — , or, more positively, a handle for 
the electrical spin manipulation4SiiL It is natural to 
expect that the presence of the spin-orbit interaction 
will influence the exchange Hamiltonian."*^ The result- 
ing corrections to the rotationally symmetric exchange 
Hamiltonian are referred to as the anisotropic exchange 
(we do not consider other sources than the spin-orbit 
interactio n^'^'^^" — ) . Stringent requirements of the quan- 
tum computation algorithms motivate studies of the con- 
sequences of the anisotropic exchange of a general form 
on quantum gates4^"— Usually, the anisotropic exchange 
is viewed as a nuisance to be minimized. ^^^^^ On the 
other hand, it was considered as a possible way of imple- 
menting the quantum gates. ^^'^^ In both views, it is of 
utter importance to know the strength and the form of 
the anisotropic exchange. Since the spin-orbit interaction 
is weak, it is enough to answer the following question: 
What is the anisotropic exchange in the leading order? 

Surprisingly, arriving at the answer was not straight- 
forward at all. The Dzyaloshinskii-Moriya^^^ interac- 
tion is of the first order in spin-orbit coupling. How- 
ever, since it couples only states split by the isotropic 
exchange, it is necessary to consider also the sec- 
ond order anisotropic exchange terms to arrive at cor- 
rect energies ^-■^ RefHi suggested such a Hamiltonian, 
which was unitarily equivalent to the isotropic exchange 
Hamiltonian, with the exchange energy renormalized in 
the second order. This was later revisited^^i^ with the 
following conclusion: In zero magnetic field, the two 
qubit Hamiltonian is, up to the second order in the linear- 
in-momcnta spin-orbit interaction, unitarily equivalent to 
the isotropic exchange Hamiltonian in the weak coupling 
limit, with the unchanged exchange energy. Further cor- 
rections appear in the third order. In the unitary op- 
erator providing the change of the basis, the spin-orbit 
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interaction appears in the linear order. These results are 
a consequence of the special form of the spin-orbit in- 
teraction, which in the leading order leads to a spatially 
dependent spin rotatioui^ 

In the short version of this article?^ we derived the 
leading order anisotropic exchange terms which appear 
in a finite magnetic field. We derived all anisotropic ex- 
change parameters in a form valid for arbitrary interdot 
coupling. We also compared the results obtained using 
the first order versus the second order treatment of the 
spin-orbit interactions. The main goal of the present 
work is a detailed assessment of the quantitative relia- 
bility of the presented anisotropic exchange model com- 
paring with exact numerical results. Specifically, we ex- 
amine the model in the strong and weak coupling regimes 
[corresponding to single (Sec. IIII|) and double (Sec. ITV)) 
dots, respectively] and in zero and finite perpendicular 
magnetic field. We also study the role of the cubic Dres- 
selhaus term (Sec. lIVPp . whose action does not corre- 
spond to a spatial texture (in the leading order) and 
could potentially become dominant over the linear terms, 
changing the picture considerably. In addition to that, 
we supply the derivations, not presented in the short ver- 
sion (Sec. Ill C|) and a detailed account of our numerical 
method (App. E]). 

The analytical pitfalls in evaluating the isotropic 
exchange are well known.— On top of that, the 
anisotropic exchange is a (very) small correction to the 
exponentially sensitive isotropic exchange, and therefore 
it is involved to extract even numerically. Our main 
conclusion here is that the presented analytical model 
is valid in all studied regimes. Quantitatively, the effec- 
tive parameters are usually within a factor of 2 from their 
counterparts derived from the numerically exact spectra. 
The main source of the discrepancy is the cubic Dressel- 
haus term. Surprisingly, in the most important regime 
for quantum dot spin qubits, namely the weak coupling, 
the Heitler-London approximation works great for the 
anisotropic exchange, even though it fails badly for the 
isotropic one. This finding justifies using simple analyti- 
cal formulas for the anisotropic exchange parameters. 



II. MODEL 

Our system is a two-dimensional electron gas con- 
fined in a [001] plane of a zinc-blende semiconductor 
heterostructure. An additional lateral potential with 
parabolic shape defines the double quantum dot. We 
work in the single band effective mass approximation. 
The two-electron Hamiltonian is a sum of the orbital part 
and the spin dependent part, 

Htot = Horb + ^ HsoS+ Hz,i = Horb + Hso + Hz, (1) 
1=1,2 



where the subscript i labels the two electrons. The orbital 
Hamiltonian is 

iforb= 51 (^« + ^0+^C- (2) 

1=1,2 

Here, Ti — h'^K.f/2m is the kinetic energy with the 
effective mass m and the kinetic momentum KKi = 
hki + eAj = —ih^i + eAi] e is the proton charge and 
Ai = Bz/2{—yi,Xi) is the vector potential of the mag- 
netic field B = {Bx, By, Bz). The potential V describes 
the quantum dot geometry 

V- = imw2^in{(r, - d)^, (r, + d)^}. (3) 

Here Iq = {h/moJo)^^^ is the confinement length, 2d mea- 
sures the distance between the two potential minima, the 
vector d defines the main dot axis with respect to the 
crystallographic axes and Eq — Hloq is the confinement 
energy. The Coulomb interaction between the two elec- 
trons is 

1 

He = -. r, (4) 

47reoer jri - 

where is the vacuum dielectric constant and is the 
dielectric constant of the material. 

The lack of the spatial inversion symmetry is accom- 
panied by the spin-orbit interaction of a general form 

Hso,i ^yVi ■ (Ti, (5) 

where the vector w is kinetic momentum dependent. In 
the semiconductor heterostructure, there are two types 
of spin-orbit interactions. The Dresselhaus spin-orbit in- 
teraction, due to the bulk inversion asymmetry of the 
zinc-blende structure, consists of two terms, one linear 
and one cubic in momenturo^^ 

wz,,, = Jc{Kl,){~Kx^^,Ky,„0), (6) 
^D3., = jc/2{Kx,,Kl^,-Ky.,Kl^,0)+Ii.c., (7) 

here H.c. denotes the Hermitian conjugate. The in- 
teraction strength 7c is a material parameter, the an- 
gular brackets in wd denote the quantum averaging in 
the z direction. Since both electrons are in the ground 
state of the perpendicular confinement, we have {K'^ i) = 
— i-^z)' value depending on the confinement 
details. A confinement asymmetry along the growth di- 
rection (here z) gives rise to the Bychkov-Rashba term?^ 

w_Bj?,i = ctBR {Ky^i, -Kx,i, 0) . (8) 

The coupling a br of the interaction is structure depen- 
dent and can be, to some extent, experimentally mod- 
ulated by the top gates potential. Equations are 
valid for a coordinate system where the x and y axes 
are chosen along [100] and [010] directions, respectively. 
Below we use the effective spin-orbit lengths defined as 
Ibr = h^/2maBR and Id = h"^ /2m^c{K^) ■ 
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The spin is coupled to the magnetic field through the 
Zeeman interaction 



Hz^i = |m-bB • (Ti = i^iB ■ (Ti 



(9) 



where g is the effective gyromagnetic factor, /is = 
eft/2TOe is the Bohr magneton (alternatively, we use a 
renormalized magnetic moment /i) and cr is the vector of 
the Pauli matrices. 

In lateral quantum dots the Coulomb energy Ec is 
comparable to the confinement energy and the cor- 
relation between the electrons strongly influence the 
stateS i^^'^^ One can compare the energies considering 



'0 



lc{r)' 



(10) 



where the Coulomb length Ic = e^m/47reoer^^ is a mate- 
rial parameter and (r) is the mean distance between the 
electrons. In GaAs Ic ~ 10 nm, while a typical lateral 
dot has ^0 ~ 30 nm, corresponding to i?o ~ 1 meV. The 
mean length (r) is of the order of the confinement length, 
if the two electrons are on the same dot, and of the in- 
terdot distance, if the electrons are on different dots. In 
the first case, the Coulomb energy is typically 3 meV. 
In the second case (one electron per dot) the Coulomb 
interaction is typically at least 1 meV. 

The strength of the Coulomb interaction precludes the 
use of perturbative methods. Therefore, to diagonalize 
the two electron Hamiltonian Eq. ([T|), we use the ex- 
act numerical treatment, the Configuration interaction 
method. Details are given in App. Below we consis- 
tently use the notation of $ for spinor and ^ for orbital 
wavefunctions. They fulfill the equations Htot^ = 
and Horb^ = E'^ , respectively. 

We use the GaAs realistic parameters: m = 0.067TOe 
(TOe is the free electron mass), g = —0.44, = 12.9 and 
7c = 27.5 eVA . The coupling of the linear Dresselhaus 
term is ^c{Kl) = 4.5 meVA and of the Bychkov-Rashba 
term \s Ubr = 3.3 meVA, corresponding to the effec- 
tive spin-orbit lengths Id — 1.26/im and Ibr = 1.72/im, 
according to the recent experiments i^i^ We use the 
confinement energy tvjjQ = 1.1 meV, which corresponds 
to the confining length /q — 32 nm, in line with an 
experiment.®^ 



A. Unitarily transformed Hamiltonian 



spin-orbit interactions are small perturbations, compara- 
ble in strength to the Zeeman term at B = 1 Tesla. 

We consider the perturbative solution of the Hamilto- 
nian Eq. ([T}. We transform the Hamiltonian to gauge 
out the linear spin-orbit terms^^i^ (we neglect the cubic 
Dresselhaus term in the analytical models) 



Htot — > U HtotU^ — Horb + Efz + Hs 



using the operator 



J7 = exp ( -^ni • CTi - • era 



where 



^br ^br 



Id 



.0 



(11) 



(12) 



(13) 



Keeping only terms up to the second order in the spin- 
orbit and Zeeman couplings, we get the following effective 
spin-orbit interactions Hgo = El so + , where 



HfJ = J2 (-^+ + K-L,^,a,Jh) , (14) 
-(a^B X m) ■ CTi. (15) 



i=l,2 
r(2) 



1=1.2 



Here, L^^i/h = XiKy^^ - ViK^^i, and 



4mlj 



± 



(16) 



Equation (jl5p describes the mixing between the Zeeman 
and spin-orbit interactions, which is linear in the spin- 
orbit couplings. It disappears in zero magnetic field, 
where only the terms in Eq. (jl4p survive - a sum of an 
overall constant shift of 2A% and the spin-angular mo- 
mentum operators. Both of these are quadratic in the 
spin-orbit couplings. 

The point of the transformation, which changes the 
form of the spin-orbit interactions, is that the trans- 
formed interactions are much weaker (being the second, 
instead of the first order in the spin-orbit/Zeeman cou- 
plings). Of course, both Hamiltonians are equivalent, 
giving the same exact energies. However, a perturba- 
tive expansion of the transformed Hamiltonian converges 
much faster. 



Analytically, we will analyze the role of the spin-orbit 
interactions in the two-electron spectrum using the per- 
turbation theory. This approach is appropriate since the 
spin-orbit energy corrections are small compared to the 
typical confinement energy. For a GaAs quantum dot the 
ratio between the confinement length and the spin-orbit 
length lo/lso ~ 10~^-^10~^. Furthermore, for a magnetic 
field of 1 Tesla, the ratio between the Zeeman energy and 
the confinement energy is ^iB/Eq ^ 10~^. Therefore the 



B. Orbital functions symmetry 

The symmetry of the two-electron wavefunctions 5" has 
important consequences, for example, in the form of se- 
lection rules for the couplings between the states due to 
the spin-orbit interactions. The choice of the potential 
in Eq. ([3]) is motivated by the fact that for small {d 0) 
and large (d — oo) interdot distance the eigenstates of 
the single particle Hamiltonian converge to the single dot 
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solutions centered at d = and x = ±2d, respectively. 
For zero magnetic field, since the double dot potential 
does not have the rotational symmetry around the z axis, 
the inversions of the coordinate along axes of the confine- 
ment potential (x and y) are the symmetries involved. 
Indeed, the orbital Hamiltonian Eq. ([2]) commutes with 
the inversion operator and ly, [Horb, Ix.y] — 0. Fur- 
thermore [Horbjl] = 0, where / = Ixly is the inversion 
of both axes simultaneously. All these operations belong 
to the C'2v group. Accordingly, the wavefunctions trans- 
form as the functions 1, x, xy, and y, which represent 
this group. If a perpendicular magnetic field is applied, 
only the total inversion operation, / = Ixly, commutes 
with the Hamiltonian and the wavefunction is symmet- 
ric or antisymmetric with respect to the total inversion 
- this is due to the lack of Ix and ly symmetry of the ki- 
netic energy operator. The Slater determinants (the two- 
electron basis that we use in the diagonalization proce- 
dure - see App. [A| have also definite symmetries, if they 
are built from single particle states of definite symmetry 
(see App. |B]). 

We define the functions 'i>± to be the lowest eigenstates 
of the orbital part of the Hamiltonian, Horb'^± = E±'^± 
with the following symmetry, 

P«'± = ±«'±, (17) 

where P/152 = /2ffi is the particle exchange operator. 
We observe that ^± have, in addition to the particle 
exchange symmetry, also a definite spatial symmetry. In 
further we assume they fulfill 

/i/2*± = ±*±. (18) 

We point out that while Eq. (fT7| is a definition, Eq. (fT8|) 
is an assumption based on an observation. In zero mag- 
netic field Iil2'^+ = +^+, follows from the Mattis- 
Lieb theoremji For the validity of Eq. (flSl) we resort to 
numerics — we saw it to hold in all cases we studied. 

Figure [1] shows the calculated double dot spectrum at 
zero magnetic field without the spin-orbit interactions. 
The two lowest states are split by the exchange en- 
ergy J. In the single dot case {d = 0), the ground state 
is non degenerate, while the first exited state is doubly 
degenerate. Increasing the interdot distance, this degen- 
eracy is removed, as the two states have different spatial 
symmetry (x and y). The energy of the states is 
separated from the higher states by an energy gap A. 
This gap allows us to consider only the two lowest or- 
bital states when studying the spin-orbit influence on the 
lowest part of the two-electron spectrum. Indeed, in the 
double dot A is of the order of 1 meV, while the spin- 
orbit interactions are two orders of magnitude smaller. 
In the case of A = 0, the two orbital states approxima- 
tion can be improved including more states (although we 
show below this is not in fact necessary for a qubit pair 
in a circular dot). 

Without the spin-orbit interactions, the eigenstates of 
the Hamiltonian Eq. ([l} are separable in the spin and 




interdot distance [nm] / tunneling energy [|J,eV] 



FIG. 1: Two-electron energy spectrum of a double dot at zero 
magnetic field as a function of the interdot distance and the 
tunneling energy. The spatial symmetries of wavefunctions, 1, 
X, xy, and y are denoted as solid, dashed, dotted-dashed and 
dotted line, respectively. The two lowest energies are labeled; 
they are split by the isotropic exchange energy J. The energy 
separation between the lowest states and the higher exited 
states is denoted by A. 



orbital degrees of freedom. We get the four lowest states 
by supplementing ^± with spinors, forming the singlet 
and triplets: 

..,4 = {"f+S, *_To, *_T+, ^-.T.}. (19) 

Here = 1 / \/2 ( 1 1|) - ||t) ) ) is a singlet spinor built out of 
two spin-1/2 spinors. To = l/\/2(|n) + lit)), T+ = |n>, 
T_ = 14,4,) are the three possible triplets; the quantization 
axis is chosen along the magnetic field. 

The symmetry leads to selection rules for the matrix 
elements between two electron states. In zero perpen- 
dicular magnetic field, because the Lz operator trans- 
forms as xy, the singlet and triplets are not coupled, 
up to the second order in the spin-orbit interactions, 
(<i>i|i/so|$2,3,4) = 0. The only contribution is due to 
the constant K^. For non zero perpendicular magnetic 
field, the singlet and a triplet are coupled only if their 
orbital parts have the opposite spatial symmetry, due to 
the term in Eq. (|15l) . The non- vanishing matrix elements 
are listed in Table ID 
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zero perpendicular field 


finite perpendicular field 




never 


ii = h 


ni 


Jl / J2 





TABLE I: Conditions on the orbital symmetries for the ma- 
trix elements (*i|0|'I'2) to be non-zero. The orbital symme- 
tries are defined by /\['i,2 ~ ji,2*I'i,2. 
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C. Effective Hamiltonians 

Here we derive effective four level Hamiltonians, which 
provide understanding for the numerical results. We fol- 
low two different approaches: (i) restriction of the to- 
tal Hamihonian, Eq. ([Ij, to the basis in Eq. p^ : (ii) 
including higher excited states through a sum rule us- 
ing the Schrieffer- Wolff transformation with the unitary 
operator, Eq. (jl2p . Then we compare the two models, 
including their simplifications using the Heitler-London 
approximation, to demonstrate the quality of their de- 
scription of the two-qubit subspace. 

We restrict the Hilbert space of the double dot to the 
four lowest functions Eq. (jl9p to describe the qubit pair. 
We start with the case of zero spin-orbit interactions. 
In the external magnetic field, the two triplets T_|_ and 
T_ are split by twice the Zeeman energy Ez — 
The restriction of the Hamiltonian Eq. ^ to the basis 
Eq. (fT9)) produces a diagonal matrix 

H,so - d\&g{E+,E^,E^+Ez,E^ - Ez). (20) 

The standard notation is to refer only to the spinor part 
of the basis states. The matrix Eq. (pO)) can be rewrit- 
ten in a more compact way using the basis of the six- 
teen sigma matrices, {cTa^iUj3^2\a,i)=Q,x,y,z (index de- 
notes a unit matrix; for explicit expressions see App. iPj) . 
The result is the so-called isotropic exchange Hamilto- 
nian (where the constant — J/4 was subtracted) 

ff^so = (J/4)cri-cr2+MB-(cri+cr2), (21) 

where the singlet and triplets are separated by the 
isotropic exchange energy J ~ E- — E+, the only pa- 
rameter of the model. 

The Hamiltonian Eq. (PT|) describes the coupling of the 
spins in the Heisenberg form. With this form, the SWAP 
gate can be performed as the time evolution of the sys- 
tem, assuming the exchange coupling J is controllable. 
The isotropic exchange has already been studied analyt- 
ically, in the Heitler-London, Hund-MuUiken, Hubbard, 
variational and other approximations, as well as numeri- 
cally using the finite-difference method. Usually analyti- 
cal methods provide a result valid within certain regime 
of the external parameters only and a numerical calcula- 
tion is needed to assess the quality of various analytical 
models. 

When the spin-orbit interactions are included, ad- 
ditional terms in the effective Hamiltonian appear, as 
the matrix elements due to the spin-orbit interactions 
{Haniso)ij {^i\Hso\^j)- Sclcction rulcs in Tab. H] re- 
strict the non-zero matrix elements to those between a 
singlet and a triplet, 

/ 2w,^ -V2u V2v \ 
, ^ 2w: 

-V2u* 0'^' 
\ V2v* / 



Here u = {wx + ^Wy), v = {wx — iwj^) and 

w = (*+|wi|*_), (23) 

where vector w is defined by the spin-orbit interactions 
Eq. ([S]). Using the sigma matrix notation, Eq. can 
be written as (see App. [DJ 

Kmso = a' • (,Ti - era) + b' • (cti x (Ts), (24) 

where the a' and b' are the spin-orbit vectors defined as 

a' = Re(*+|wi|1'_), (25a) 
b' = Im(*+|wi|'J'_). (25b) 

The standard exchange Hamiltonian follows as 

^ex ~ Hiso + Hg^^^g^, (26) 

and we refer to it in further as the first order (effective 
model) to point the order in which the spin-orbit inter- 
actions appear in the matrix elements. Note that we 
repeated the derivation of ReflssI additionally including 
the external magnetic field. As we will see below, com- 
parison with numerics shows that treating the spin-orbit 
interactions to the linear order only is insufficient. 

To remedy, we generalize the procedure of Ref. i59| to 
finite magnetic fields. This amounts to repeating the 
derivation that lead to Eq. 1201 this time starting with 
the unitarily transformed Hamiltonian Eq. (jlip . In this 
way, the linear spin-orbit terms are gauged out and the 
resulting effective Hamiltonian treats the spin-orbit inter- 
actions in the second order in small quantities (the spin- 
orbit and the Zeeman couplings). The transformation as- 
serts that the original Schrodinger equation Htot^ = E^ 
can be equivalently solved in terms of the transformed 
quantities Htot{U^) = E{U^), with the Hamiltonian 
H = UHtotU''. The transformed Hamiltonian H is the 
same as the original, Eq. ([I}, except for the linear spin- 
orbit interactions, appearing in an effective form Hgo- 
We again restrict the basis to the lowest four states and 
for the spin-orbit contributions we get 

($,|iJ,o|$,). (27) 

Using the selection rules and the algebra of the Pauli 
matrices, we get the exchange Hamiltonian (for obvious 
reasons, we refer to it as the second order model) 

Hex = ( J/4)cri • (72 + fi{B + Bso) • (o-i + (72) 

+ a • (cTi - (72) + h- {(7i X 0-2) - 2K+. 

Compared to the first order model Eq. ([M)) . the func- 
tional form of the second order model Hamiltonian is the 
same, except for the effective spin-orbit magnetic field 

fxBso = z(A'_/?i)(*_|i,,i|*_), (29) 

which appears due to an inversion symmetric part of iJ , 
Eq. (1141) . The spin-orbit vectors, however, are qualita- 
tively different 

a = X Re(*+|ni|^'_), (30a) 
b = X Im(*+|ni|*_). (30b) 
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We remind that the second order effective model Hamil- 
tonian Eq. (1^ refers to the four fmictions in Eq. 
unitarily transformed {[/$i}i=i,...,4. The agreement be- 
tween the second order effective model and the numerical 
data is very good, as we will see below. 



D. First order effective Hamiltonian in zero field 

In this section we give H'^^ explicitly for zero B and 
diagonalize it. This is the only case for which is possible 
to give an analytical solution. For zero magnetic field, 
one can choose the functions to be real. Then the 
matrix elements of the spin-orbit operator w in Eq. ([5]) 
are purely imaginary and a' = . With the spin quan- 
tization axis chosen along the vector b', the 4x4 matrix, 
Eq. (pS)) . takes the form of 



TT' _ 

ex 



I -ZJI^ 2i6' \ 
-2i6' J/4 
J/4 

\ J/4 y 



(31) 



The upper left 2x2 block of this matrix is a Hamilto- 
nian of a spin 1/2 particle in a fictitious magnetic field 
S = (0, 2&', J/2)//i. The eigenstates of this Hamiltonian 
are spins oriented along the magnetic field B. Since the 
matrix in Eq. ((3T|) is block diagonal, it is easy to see it 
can be diagonalized with the help of the following matrix 



/O 1 o\ 
10 


VO 0/ 



(32) 



The Hamiltonian Eq. (|3T|) can be diagonalized by 



Hdiag 



f ~J/4~\fj,B\ \ 

-J/4:+\fiB\ 

J/4 

V J/4:J 



, (33) 



where = 4(6')^ -I- j2/4. In the notation of the Pauh 

matrices, (see App. ID|) . 

e = exp (^^i^^ = (^-^^i'^B.i " <^B.2)^ , (34) 

where tan 6* = Ab' / J and as = cr ■ B/B. 

The unitary transformation in Eq. (j34]) performs 
the rotation of the two spins in the opposite sense. The 
Hamiltonian can be interpreted as a rotation of the elec- 
tron around a spin-orbit field when transferred from one 
dot to the otherj^ The spectrum given by Eq. ([33)) qual- 
itatively differs from the numerics, which shows there is 
no infiuence on the exchange in the second order of the 
spin-orbit couplings. 




0.4 0.6 0.8 

magnetic field [T] 



1.0 



FIG. 2: Two-electron energy spectrum of a single dot in per- 
pendicular magnetic field. The lowest states are labeled by 
the total angular momentum L. The Zeeman and spin-orbit 
interactions are neglected. The two regions marked by boxes 
are magnified on Fig. [3] and Fig. [l] 



III. SINGLE DOT 

We start with the single dot case, corresponding in 
our model to d = 0. The analytical solution of the single 
particle Hamiltonian T + V is known as the Fock-Darwin 
spectrum. The corresponding wave functions ip and the 
energies e are 

Vn/(r„^,)=Cpl'le-''^/2Li;i(p,2)e''^% (35) 



{2n+\l\ + l) + B^l, (36) 
Zra 



where p, = rjls and Ib = [l^'^ + {eB^/2hf]-^/'^ is the 
magnetic length; n and / are the radial and the angular 
quantum numbers, C is the normalization constant and 
L^n are the associated Laguerre polynomials. 

Let us consider now the orbital two electron states 
^, eigenstates of iJorb, Eq. 0. The Coulomb operator 
He commutes with the rotation of both electrons around 
the z axis, that is, the Coulomb interaction couples only 
states with the same total angular momentum. This al- 
lows us to label the states with the quantum number 
L — L1+L2, the total angular momentum. Furthermore, 
the Hamiltonian -fforb commutes with any spin rotation 
of any of the electrons, which expresses the fact that the 
Coulomb interaction conserves spin. Therefore we can 
consider the full two electron wavefunctions obtained by 
supplementing the orbital part ^ with a spinor, respect- 
ing the overall wavefunction symmetry, similarly as in 
Eq. dUl). 

The two-electron spectrum, without the Zeeman and 
the spin-orbit interactions, is shown in Fig. [21 At zero 
magnetic field the ground state is a non-degenerate sin- 
glet state with total angular momentum zero L — 0. The 
next two degenerate states are triplets with L — ±1 and 
their degeneracy is split by the magnetic field. Focus- 
ing on the two lowest states, most relevant for the qubit 
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pair, they cross at B « 0.43 T, so one can turn the 
ground state from the singlet to the triplet by applying 
an external magnetic field. In the presence of spin-orbit 
interactions, the crossing is turned into anticrossing, as 
described below. 



A. Spin-orbit correction to the energy spectrum in 
magnetic field 

Suppose some parameter, such as the magnetic field, 
is being changed. It may happen at some point that 
the states of the opposite spin become degenerate. Such 
points are called spin hot spots. Here, because of the 
degeneracy, weak spin-orbit interactions have strong ef- 
fects. For the spin relaxation, spin hot spots play often 
a dominant rolci^ 

We are interested in the changes to the spectrum due 
to the spin-orbit interactions. Let us neglect the cubic 
Dresselhaus in this section. To understand the spin-orbit 
influence, it is important to note the following commuta- 
tion relations for the linear spin-orbit terms 



[W£)4 • CTi + wd,2 ■ cr2, J-] = 0, 



(37) 



where J± = '^i{Lz,i ± Sza)- These commutation rules 
hold for any magnetic field B. Since the Hamiltonian 
Eq. ^ commutes with the operator J±, we can label 
the states using the quantum numbers J+ = L -t- 5"^ and 
J_ = L — Sz- The spin-orbit interactions couple only the 
states with the same quantum numbers J-|_ and J_, for 
Bychkov-Rashba and Dresselhaus term, respectively. 

Let us focus on the part of the spectrum close to i? = 
and on the states with L = ±1, Fig. [21 The degeneracy 
of the states is removed by the spin-orbit interactions, as 
shown in Fig. [S] 

Let us now use the Hamiltonian Eq. (|lip . to under- 
stand the influence of the spin-orbit interactions. The 
degeneracy of the states with angular momenta L = ±1 
makes the description with the lowest two orbital states 
questionable. Therefore now we take 3 orbital states and 
repeat the derivation of the second order effective Hamil- 
tonian, obtaining a 7 x 7 matrix. The basis functions 
are 



{^^}^=l,...,7 = + *-ro, *_T+, *_r_, 

^'_To,'i''_T+,^'_T^}, 



(38) 



where is the electron wavefunction with angular mo- 
mentum L — 0, and and have angular momentum 
L — +1, and L = —1, respectively. Since the magnetic 
field is negligible with respect to the spin-orbit couplings, 
the Hamiltonian Eq. ([T5|) is negligible. Because of the se- 
lection rules. Tab. U the contributions from Eq. (|T4l) in 
the basis Eq. psp . gives non zero matrix elements only 
for the following pairs, -T±\H so\'^ -T±) = ±K_, and 
(^''_r±|i?,o|*'_T'±) = ±K_. For the GaAs parameters. 
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FIG. 3: (Color online) Magnified region from Fig. O En- 
ergy spectrum of a single dot for small perpendicular mag- 
netic field. Only the states with the total angular momentum 
L — itl are plotted. A constant shift is removed from the 
spectrum. Each state is labeled by the quantum numbers 

{J+,J-,T^). 
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FIG. 4: (Color online) Lowest energy levels in the anticrossing 
region marked in Fig. [2] A constant shift was removed from 
the spectrum. The quantum numbers (J+,J_,Ei) label the 
states. Insets show the anticrossing regions. 



— 0.16/1 eV. In the region of small magnetic field, the 
states with J± — are coupled by the spin-orbit inter- 
actions and the lifting is in the second order in the spin- 
orbit couplings. The other states are not coupled since 
they have different values of J±. Therefore we conclude 
that the 2-orbital state approximation can be used also 
for the single dot case (or strongly coupled double dots) , 
because the spin-orbit interactions do not mix the states 
'i' and 4'' in the basis Eq. (|38l) . Note that as the coupling 
is forbidden by the inversion symmetry, the claim holds 
for an arbitrary oriented magnetic field. 

Let us now discuss the second degeneracy region 
marked in Fig. [21 magnified in Fig. [31 The spin-orbit 
interactions induce two anticrossings. The first is due 
to the Bychkov-Rashba term, since the crossing states 
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have different J_, but the same J+ = and couples 
the singlet S and triplet T_(_. The second is due to the 
Dresselhaus term which couples states with J_ = 0, the 
singlet S and the triplet T_ . The central point is a cross- 
ing point, because the crossing state differ in both J+ 
and J- . The splitting energy can be evaluated using the 
unitarily transformed Hamiltonian Eq. (jlip . Using the 
degenerate perturbation theory, one can estimate ana- 
lytically, using Eq. (fMl) and Eq. (|T5|) . the value of the 
two gaps to be Abr ~ 4:^/2^BIo/Ibr = O.lbfxeV and 
Ad ~ 4:^/2fiBlo/lD — 0-58fieV. These values are consis- 
tent with the numerical values. 



and the spin-orbit induced magnetic field 

^-WlU'rW^ (42) 

-{^L\LAi^R){i^RHL)). 

The explicit formulas for the vectors in Eqs. (|40|) - (j42|) 
are in App. [Cl Differently from the spin-orbit vectors 
in Eq. (HUl) . the vectors in Eq. (HIT) reveal explicitly the 
anisotropy with respect to the magnetic field and dot 
orientation^SiS, (note that x and y in the definition of n, 
Eq. (fT3|) are the crystallographic coordinates). 



IV. DOUBLE DOT 

The double dot denotes the case when the interdot 
distance is of the order of the confinement length. In the 
next sections we discuss our effective models, Eq. (pS)) 
and Eq. (|26p in the double dot regime and compare them 
with numerics. 



A. Heitler-London approximation 

The analytical solution for the two electron wavefunc- 
tions in a double dot potential is not known. We consider 
here the Heitler-London ansatz since it is a good approx- 
imation at large interdot distances and we can work out 
the spin-orbit influence on the spectrum analytically. For 
this purpose, we compute the spin-orbit vectors, Eq. ((25)) 
and Eq. (ISH]) . for our models. 

In the Heitler-London ansatz, the two electron eigen- 
functions are given by 

1 



v/2(l±|(^L,i|V'fl,i)P) 



*± = ,,,, ,, ,,.,, (IV'£,l)IV'fl.2)±|V'fl.l)IV'L.2)), 

(39) 

where | V'L(i?),i) is a single electron Fock-Darwin state cen- 
tered in the left (right) dot occupied by the z-th electron. 
Below, in Eqs. (HHjl - lH^ . we skip the particle subscript i, 
as the expressions contain only single particle matrix ele- 
ments (all ip, w, n, Lz would have the same subscript, say 
i — \). With this ansatz, the spin-orbit vectors, Eq. 
follow as 



b' 



^1-\{^l\^rW 
i 



(?Al|w|i/>l), 



(40a) 



:(VL|w|7/;fl)(V'fl,|V'L). (40b) 



Similarly we get the spin-orbit vectors, Eq. ([30|) . as 



a = 



= (?Al|B X n|VL), (4l£ 



:(VL|Bxn|Vi?)(V'i?|V'L>(41b) 



B. Spin-orbit correction to the energy spectrum in 
zero magnetic field 

In the previous sections, we have derived two effec- 
tive Hamiltonians, H'^^, and H^x, given by Eqs. 
and Eqs. ((24|) - p6)) . respectively. We now compare the en- 
ergy spectrum given by these models with exact numer- 
ics. We present the spin-orbit induced energy shift, the 
difference between a state energy if the spin-orbit interac- 
tions are considered and artificially set to zero. For each 
model we examine also its Heitler-London approxima- 
tion, which yields analytical expressions for the spin-orbit 
vectors, as well as the isotropic exchange energy (given in 
Sec. IV. A and Appendix C). Thus, the effective models 
in the Heitler-London approximation (we denote them by 
superscript HL) are fully analytic. The non-simplified 
models (we refer to them as "numerical") require the 
two lowest exact double dot two-electron wavefunctions, 
which we take as numerical eigenstates of Horb- 

Apart from the energies, we compare also the spin- 
orbit vectors. Since they are defined up to the relative 
phase of states '5+ and ^f-, the observable quantity is 
c' = y^{a'Y + {b'Y a-iid analogously for c = ^a? + P. We 
refer to these quantities as the anisotropic part of the 
exchange coupling. 

Figure O shows the spin-orbit induced energy shift as 
a function of the interdot distance for each of the four 
states. The exact numerics gives a constant and equal 
shift for all 4 spin states, with value — 0.54/zeV. Let us 
consider the second order model, Eq. (pS)) . For zero mag- 
netic field, all spin-orbit vectors are zero, as is the effec- 
tive magnetic field. The only contribution comes from 
the constant term = — 0.54/ieV that is the same for 
all states. Our derived spin- model, Eq. (pS)) . accurately 
predicts the spin-orbit contributions to the energy. On 
the other hand, the first order models H'J^^ and H'^^ 
are completely off on the scale of the spin-orbit contribu- 
tions. The exchange Hamiltonian H'^^ does not predict 
the realistic spin-orbit influence on the spectrum, even in 
the simple case when the magnetic field is zero. 

Figure [6] shows the non zero parameters for all four 
models. The exact isotropic exchange J decays exponen- 
tially with the interdot distance. The same behavior is 
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FIG. 5: Spin-orbit induced energy shifts at zero magnetic field 
as a function of the interdot distance. Exact numerics (solid), 
first order model H'^^ (dotted) and first order model in HL 
approximation (dashed) are given, a) Singlet, b) triplet To, 
c) triplet T+, d) triplet T-. The results of the second order 
model (both H^J" and Hex give the same) are indiscernible 
from exact numerical data. 







1 


1 ' 1 ' 1 ' . 


10-' 




^ — 




10-^ 




c' ""^^^ 


N. \ 


10-= 






\\ ' 










10-^ 




1 


* \\ 
\\ 

■V \\ 

^ \ \ 

1 , 1 , 1 \,\\ ■ 



20 40 60 80 



interdot distance [nm] 
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London approximation for the isotropic exchange (solid) and 
the anisotropic exchange of the first order model (dashed). 



predicted in the Heitler-London approximation. It de- 
cays exponentially, but deviates from the numerical re- 
sults. As for the anisotropic exchange, the first order 
model H'^^ gives an exponentially falling spin-orbit pa- 
rameter c', an order of magnitude smaller than J. In 
contrast, the second order model H^x predicts zero spin- 
orbit anisotropic exchange. First main result, proved nu- 
merically and justified analytically by the Hamiltonian 
Hex^ is that at zero magnetic field the spin-orbit vectors 
vanish^ up to the second order in spin-orbit couplings at 
any interdot distance. In the transformed basis, there is 
no anisotropic exchange at the zero magnetic field due the 
spin-orbit interactions, an important result for the quan- 
tum computation. Indeed, since the exchange energy can 



FIG. 7: Spin-orbit induced energy shifts at 1 Tesla perpen- 
dicular magnetic field versus the interdot distance, a) Energy 
shift of the Singlet 5* in the exact numerics (solid) is compared 
to the numerical (dashed) and the Heitler-London approxima- 
tion (dotted) first order model. In b) similar comparison is 
made for the second order model. Panels c-d) are analog of 
a-b) showing the energy shifts of the triplet T+. 

be used to perform a SWAP operation, this means that 
the spin-orbit interactions do not induce any significant 
errors on the gate operation. The only difference is the 
computational basis, which is unitarily transformed with 
respect to the usual singlet-triplet basis. 

C. Finite magnetic field 

In the presence of a perpendicular magnetic field the 
structure of the spin-orbit contributions are quite dif- 
ferent with respect to the zero field case. First of all, 
anticrossing points appear, where the energy shift is en- 
hanced. Figure [7] shows the spin-orbit contributions in 
a finite magnetic field . We plot only the anticrossing 
states, the singlet S and the triplet T+. The prediction 
of the first order model is shown in the left panels of 
Fig. [T] As in the case of zero magnetic field, this model 
is off from the numerical results. In particular, it still pre- 
dicts a zero contribtition, except close to the anticrossing 
point. We note that the discrepancy is not connected to 
(a failure of) the Heitler-London approximation, as using 
the exact numerical two electron wavefunctions does not 
improve the model predictions. 

In the right panels of Fig. [71 the comparison between 
the second order model and the numerics is provided. 
We observe that the model is very close to the numerics, 
even though the Heitler-London approximation predicts 
the crossing point in a different position. The predictions 
of the numerical second order model H^x is consistent 
with the exact numerics. The only discrepancy is due to 
the infiuence of the cubic Dresselhaus term, as we will 
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FIG. 9: Spin-orbit induced energy shifts of a double dot sys- 
tem with interdot distance of 55 nm versus the perpendicular 
magnetic field, a) singlet S, b) triplet To, c) triplet T+, d) 
triplet T-. Exact numerics (solid) and the numerical second 
order model Hex (dashed). 



see in the next section. 

To get more insight, in Fig. [S] we have plotted the pa- 
rameters of the models. Fig. [8^ shows the anisotropic 
exchange strengths in the two models. The first order 
model H'^^ predicts the anisotropic exchange decreasing 
with the interdot distance, similar to the isotropic ex- 
change energy. For large interdot distance the anisotropic 
exchange c' disappears. This means there is no influence 
on the energy due to the spin-orbit interactions. On the 
other hand, for the second order model Hex the conclu- 
sion is different. For large interdot distances c^^ and c 
are linear in d. Furthermore, the anisotropic exchange 
computed in the Heitler-London ansatz is very close to 
the numerical one. We make a very important obser- 
vation here: surprisingly, concerning the anisotropic ex- 
change the Heitler-London is quite a good approxima- 
tion for all interdot distances even in a finite magnetic 
field. Therefore, despite its known deficiencies to evalu- 
ate the isotropic exchange J, it grasps the anisotropic ex- 
change even quantitatively, rendering the spin-orbit part 
of the second order effective Hamiltonian Hex fully ana- 
lytically. One can understand this looking at Eqs. (pO|. 
The anisotropic exchange vectors are given by the dipole 
moment of the matrix element between the left and right 
localized state (see App. C for explicit formula). This 
dipole moment is predominantly given by the two local 
maxima of the charge distribution (the two dots) and is 
not sensitive to the interdot barrier details, nor on the 
approximation used to estimate the lowest two orbital 
two-electron states. This is in strong contrast to the 
isotropic exchange, which, due to its exponential char- 
acter, depends crucially on the interdot barrier and the 
used approximation. 



Figure [8b shows the isotropic exchange J, and the ef- 
fective magnetic field induced by the spin-orbit interac- 
tions iiBso compared to the Zeeman energy 2/ii?. We 
see the failure of the Heitler-London approximation for 
J. Although the numerical calculation and the analyt- 
ical prediction have the same sign (this means that the 
ground state is the triplet) they differ by an order of 
magnitude. The Zeeman energy is constant and always 
much larger than the effective spin-orbit induced mag- 
netic field fiBso- Consequently, the effective field can be 
always neglected. The point where the Zeeman energy 
equals to the isotropic exchange (close to d = 18nm) is 
the anticrossing point, where the spin-orbit contributions 
are strongly enhanced, as one can see in Fig. [71 

Let us consider a double dot system at fixed interdot 
distance of 55nm, corresponding to zero field isotropic 
exchange of lfj,eV. In Fig. [Hi the spin-orbit contributions 
versus the magnetic field are plotted for the second order 
model Hex and the exact numerics. We can conclude 
that to describe the spin-orbit infiuence on the states in 
a double-dot system it is important to use the second 
order Hamiltonian Hex- 

In Fig. [TUj the spin-orbit parameters versus the mag- 
netic field are plotted. The main infiuence on the spin is 
due to the Zeeman interaction in the whole range of B, 
since /ii?so is several orders of magnitude smaller than the 
Zeeman energy. At the ground state anticrossing point, 
the isotropic exchange crosses zero, while the anisotropic 
parameter c is finite, leading to spin hot spots. Apart 
from these, since the anisotropic exchange is two orders 
of magnitude smaller than the Zeeman energy, the spin- 
orbit induced energy shifts are minute. 
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FIG. 11: The spin-orbit induced energy shift as a function of 
the interdot distance (left panels) and perpendicular magnetic 
field (right panels), a) Singlet in zero magnetic field, c) singlet 
at 1 Tesla field, b) and d) singlet and triplet T+ at 55 nm. 
The numerical second order model H^x (dashed line), exact 
numerics (dot-dashed line) and exact numerics without the 
cubic Dresselhaus term (solid line) . 



D. Cubic Dresselhaus contributions 



V. CONCLUSIONS 

We analyzed the spin-orbit influence on two electrons 
confined in a lateral double quantum dot. We focused on 
the lowest part of the Hilbert space, which corresponds 
to a qubit pair. In Ref. l62l a Hamiltonian for such pair 
was proposed, with the spin-orbit interactions giving rise 
to an anisotropic exchange interaction. Within a uni- 
tarily transformed basis, this interaction is encoded into 
two real three dimensional spin-orbit vectors. These, to- 
gether with the isotropic exchange energy and the mag- 
netic field vector, completely parametrize an effective two 
qubit Hamiltonian. In this work, we examined the quan- 
titative validity of this effective Hamiltonian. 

In addition to a numerical study, we also provided 
the details of the effective Hamiltonian derivation, which 
were skipped in Ref. l62!. We noted that it can be di- 
agonalized analytically if the effective spin-orbit vectors 
are all aligned with the external magnetic field, the only 
exactly solvable case (apart from the trivial case of no 
spin-orbit interactions present). We also evaluated the 
spin-orbit vectors in the Heitler-London approximation 
and compared the analytical results with their exact nu- 
merical counterparts. 

There are three possible sources for a discrepancy be- 
tween the model and the exact data: the higher excited 
orbital states of the quantum dot, the higher orders of 
the effective (unitary transformed) spin-orbit interactions 
and the cubic Dresselhaus term. Elucidation of their im- 
portance is one of the main results of this work, i) We 
find the cubic Dresselhaus term is the main source of the 
discrepancy. In a typical double dot regime and a mod- 
erate field of 1 Tesla, it brings an error of ~ 0.1/ieV for 
the energies, while the two other mentioned corrections 
have an order of magnitude smaller influence, ii) We 
find the effective Hamiltonian describes both the weak 
and the strong coupling regimes (the single dot repre- 
sents the strongest possible coupling), iii) Surprisingly, 
the spin-orbit vectors obtained within the Heitler-London 
approximation are faithful even at a finite magnetic field. 
Overall, we find the anisotropic exchange Hamiltonian to 
be generally reliable, providing a realistic and yet simple 
description for an interacting pair of spin qubits realized 
by two coupled quantum dots. 



Finally we consider the role of the cubic Dresselhaus 
term. The Schrieffer- Wolff transformation does not re- 
move it in the linear order. Figure [TT] shows the en- 
ergy shifts induced by the spin-orbit interactions also in 
the case where we do not take into account the cubic 
Dresselhaus term. One can see a very good agreement 
between the second order model Hf.^ and the exact nu- 
merics where the cubic Dresselhaus term was omitted. 
Therefore we can conclude that the main part of the dis- 
crepancy we see in the spin-orbit induced energy shifts 
are due to the cubic Dresselhaus term. 
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Appendix A: Numerical method 



triplets To, and T_. We denote the new states by 



Here we describe the numerical method we use to diag- 
onalize the two electron Hamiltonian Eq. ([T|). We proceed 
in three steps 7^ We first diagonalize the single electron 
Hamiltonian H = T + V, using the numerical finite dif- 
ferences method with the Dirichlet boundary condition 
(vanishing of the wave function at boundaries). We do 
not consider the spin dependent part (spin-orbit, Zee- 
man) at this step. This allows us to exploit the symme- 
tries of the confinement potential. The single electron 
Hamiltonian is diagonalized by the Lanczos method^. 
The typical number of points in the grid we use is 60 x 60, 
giving relative precision of the energy of order 10~^. 

In the second step, using the obtained single elec- 
tron eigenstates Ci)}, we construct the two electron 
states. We use them as a basis in which the two elec- 
tron orbital Hamiltonian Eq. ([2]) is diagonalized. The 
two-electron states are constructed as symmetric 

1*^)) = i=(|V',,i)l^.,2) + |^,M)|V^.,2)) for *^iM) 



(A2) 



and antisymmetric 



1 



^/2 



{\^P^.lM.2)-\^J^m,2)), (A3) 



with respect to the particle exchange. We choose ris.e. 
single electron orbitals, typically rig.e. = 21. The total 
number of the two particle states is then ri^ ^ . 

The spatial symmetry allows us to reduce the dimen- 
sion of the two electron Hamiltonian matrix to diago- 
nalize. Namely, the matrix is block diagonal, with the 
basis functions grouped according to the spatial symme- 
try (1, x,y, xy) and particle exchange symmetry (±1). 
This results in 8 blocks and holds for zero perpendicular 
magnetic field. In a finite field, we get 4 blocks, as there 
are only two spatial symmetries possible (1, and x). Each 
block is diagonalized separately. 

The matrix element of the two-electron Hamiltonian, 
Eq. ([2]), in our basis is 



+5aM J dri J dr2*i^^^) 



1 



AneoSr |ri - 



(n,m) 



(A4) 



The last term in Eq. ()A4|) is due to the Coulomb interac- 
tion and it leads to off diagonal terms in the Hamiltonian. 
We diagonalize the matrix defined in Eq. (jA4[) to get the 
eigenspectrum {(5'i,i?i)}. 

In the third step, we add the spin dependent parts to 
the Hamiltonian. We construct a new basis by expand- 
ing the wavefunctions obtained in the previous step by 
the spin. The orbital wavefunction VP^ gets the spinor 
according to its particle exchange symmetry. The sym- 
metric function gets the singlet S while the antisymmet- 
ric appears in three copies, each with one of the three 



1$, 



I*,;) IS) 



(A5) 



where |E) corresponds to one of the 4 spin states. The 
matrix elements of the total Hamiltonian Eq. ([T]) are 



2/i|B|(fe,T+ -fe,T_))'^i,^fe,S' + 



+ ^ (vI/,|w,|vI/,,)-(Sk,|I]'), 



(A6) 



where the last term is the matrix element of the spin-orbit 
interactions. The resulting matrix is diagonalized to get 
the final eigenstates. We choose a certain number Ug of 
lowest ^'i states, depending on the required precision. In 
our simulations = 250, resulting to the accuracy of 
the order of 10~^ meV for the energy. 



Coulomb integral 

Computationally most demanding are the Coulomb in- 
tegrals. Indeed, the typical size of the Hamiltonian ma- 
trix, in the second step, is 441 x 441, requiring at least 10^ 
Coulomb integrals. Writing functions involved in the Eq. 
(|A4[) as Slater determinants, we can express the integral 
as a sum of terms such as the following 



Cijki 



dYidY2 



V'»(ri)*V'i(r2)*V'fc(ri)V';(i"2) 



|ri -ral 



dYidV2- 



(A7) 



where Fik{v) — ijji{r)*t/jk{r). The symmetry of the 
Coulomb integral Cijki — Cjuk reduces the number of 
needed matrix elements to a half. For the single dot, ipi 
are the Fock-Darwin functions and it is possible to de- 
rive an analytical formula for Cijki- In our case, since the 
single particle functions are given numerically, we have 
performed a numerical integration. Using the Fourier 
transform, we can reduce the 4-dimensional integration 
to two dimensional 



Aireoer J |q| 



where 



27r 



drTik{r) cxp(iq • r). 



(A8) 



(A9) 



For the evaluation of the Fourier transforms, we use the 
Discrete Fourier Transform algorithm with the attenua- 
tion factors, as described in Refill. 
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We compute 
mula 



according to the perturbative for- 



N^.Ny ki+k2<N kiM 



Gijki — 



n,m fci,A;2=0 (1,(2=0 



/(?i,/2,n,m); 



(-n Vfci-'i) f_„ \{k2-i2) , , 



(AlO) 



where /(q)|g„„ — J'ik{ci)J^ji{~(l) is calculated in the 
point qnm, N is the perturbative order (the order of the 
Taylor expansion), and Ny are the number of grid 
points in the x and in the y direction, respectively. The 
coefficients I2, n, m) depend only on the geometry of 
the grid and are defined as 



I{li,l2,n,m) 



dx / dy 



y- 



(All) 



Here ilx = ^/'2)Sx, {n + l/2)6x) is the integration 

region and is the grid spacing along x. Similarly for 
the y direction. In our simulations we use the previous 
formula up to the 4-nd order in the Taylor expansion. 
The achieved relative precision is 10~^, with the compu- 
tational time for one Coulomb element ~ 50 ms. 



Appendix B: Two electron symmetry 

Suppose the single particle Hamiltonian commutes 
with certain set of operators {Oa}, and therefore the sin- 
gle particle states ipi can be chosen such that they have 
definite symmetries forming a representation of the group 
O of the symmetry operators 



(Bl) 



For example, since the double dot potential has inver- 
sion symmetry along x axis, is in the group O, while 
o!p = ±1 - the states are symmetric or antisymmetric 
with respect to x inversion. Now consider the two elec- 
tron states Eq. (|AHIA3I) . These states also have 
definite symmetry if a certain operator from O acts si- 
multaneously on both particles 



(B2) 



For our case of the symmetry group C2V, since — 
±1, the set of all possible products of two characters is 
the same as the set of characters for a single particle, 
{Oa(>L}i,j = {Oa}i- This mcans the two particle states 
will form the same symmetry classes as single particle 
states with the same characters. 



Appendix C: Heitler-London approximation 

In the Heitler-London approximation, the exchange en- 
ergy is calculated as 

JhL = (*-|i/or6|*-> - (*+|i/orb|*+> (CI) 



with the functions ^± given in Eq. (|39p . The single par- 
ticle ground state wavefunction of the Fock-Darwin spec- 
trum is 



V'oo(a;,2/) 



1 



exp 



y 



2P 



(C2) 



where Ib is the effective confinement length defined by 
1% = III y/T~\^We%jAW . The wavefunctions V'L(fl) 
are obtained shifting the Fock-Darwin ground state to 
(ztZod, 0). In the presence of the magnetic field we have 
to add a phase factor because of the gauge transforma- 
tion A' = B/2{-y,x±d) A = B/2{-y,x)] we have 



ipL{R) = exp 



±idC^ 



'ipooix±lod, y). 



(C3) 



The overlap between the left and right functions is 

^ = (^L I V-i?) = exp [-(d^ + (C4) 
and the exchange energy is 
hujQ 



Jhl — 



sinh[2Cd2(l -f ,92) 
eM(:dH^]h{Cd^i)^)) 



{c,v^(exp [-Cd^]I^{Cd^) + 
1 - exp hCrf']) - 



(l-Erf(dVC)) }, 



2d 

7^ 



(C5) 

where /q is the zeroth-order modified Bessel function 
of the first kind. The factor Cg is the ratio between 
the Coulomb strength and the confinement energy, Cs = 

yji: /2 /ATTeoErlo fiwo . Similar formula can be found in 
Ref.74 for a quartic confinement potential. The for- 
mula (jC5|) has been derived in RefUl ( in the original 
paper there is a trivial typo that we correct). 

The two electron energies for the states and 
are 

Eri + Ewni ± {EcE + Ewce) 



E± = 2fiwoC + 
where 



1±J72 



(C6) 



Eri 

EWri 

EcE 

EWcE 



fuvQCsVCexp [-Cd^]lo{Cd^), 



hioc) 



2d 



2d\l - Erf(dv/C)) - ^ exp [-(d^] 



hujoCs^/Cexp [-Cd^2 + ^^)]lo{Cd^^^) 



The components of the vectors a' and b' are 
a'x =0, a'y^ 0, 

Cd 



b' = 



2mQ~^jT^WTQ 



2mlbr Vl - ^0 



(1-^2). 



(C7) 
(C8) 

(C9) 
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where 4r and Id are the spin-orbit lengths for the Rashba 
and Dresselhaus respectively. 

The matrix elements of the vector n are 

and 

Appendix D: Spin matrices 

In the singlet and triplet basis, one can evaluate the 
sixteen matrices which can be formed as the direct prod- 



uct of two Pauli matrices and the identity. Here we list 
only the matrices needed for our purposes, and we re- 
group them to combinations in which they appear in the 
text. 



(C12) 



/-3 \ 
10 

10 

V 1 y 



(Dl) 



CTi — 472 



-{ 



CTi X 0-2 



<Tl + 0-2 



/ ~V2 V2\ 




/ 


-721 -721 \ 




^0 


2 


o\ 
















2 








-V2 






721 













V \/2 


) 




V 


721 


J 




lo 





0/ 


/ -^/2i \/2i\ 




/ 72 72\ 




f 





21 


o\ 



















21 








721 









72 


1 











^ -V2i 





/ 




72 / 




V 








0/ 


/ \ 




/o 


^ 




/ 





\ 


72 72 







721 -721 













72 







-721 







2 





yo 72 j 




\o 


721 J 




\0 





-2/ 



}• 



(D2) 



(D3) 
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